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Описывается алгоритм расчёта параметров электротехнической схемы замещения озонатора на основе расчёта электростатиче- 
ского поля. На основе решения интегрального уравнения Фредгольма первого рода представлен расчёт электростатического по- 
ля электродной системы неканонической формы в кусочно-однородной водо-воздушной среде. При наличии капель воды ра- 
спределение диэлектрической проницаемости становится неконтролируемым, поэтому расчёт параметров эквивалентной элек- 
тротехнической модели озонатора предлагается производить с использованием зашумлённых осциллограмм тока и напряже- 
ния на основе интеграла Дюамеля - некорректного интегрального уравнения Фредгольма. 


В системах газо- и водоочистки большой инте- 
рес представляет использование импульсного 
барьерного разряда. Эффективность очистки опре- 
деляется электрофизическими параметрами систе- 
мы и особенностями протекания физических про- 
цессов в разряде, которые зависят от формы дей- 
ствующего напряжения в цепи. Поэтому важным 
является исследование систем очистки как элек- 
трофизического устройства, определение макси- 
мальных напряжённостей электрического поля 
при фиксированной конфигурации электродов, 
электроёмкостей и сопротивлений системы. 

Количественное определение указанных харак- 
теристик позволит расширить представление о ха- 
рактере и последовательности физических процес- 
сов в разрядном промежутке и осуществить пра- 
вильный выбор оптимальных электрических пара- 


метров, что в свою очередь обеспечит необходимые 
эксплуатационные свойства системы. Таким обра- 
зом, на основе исследования электрофизических 
процессов можно осуществить оптимизацию элек- 
трической системы, снизить энергозатраты, улуч- 
шить технологию очистки воды. 

При феноменологическом описании электри- 
ческих разрядов широко используется их описание 
как объектов электрической цепи [1-4]. В основе 
таких подходов лежит замена электрофизических 
явлений, происходящих в разряде, явлениями, 
протекающими в электрической цепи, состоящей 
из известных электротехнических элементов. 

При исследовании физики разрядов наиболее 
доступными для измерения являются напряжение 
и(1) и ток і(і) в одноконтурной цепи с разрядным 
промежутком. С помощью измеренных значений 
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и(1 ) и /(/) можно определить параметры двухполю- 
сника, который является эквивалентной электро- 
технической моделью озонатора (рис. 1, й). Здесь 
С 6 - ёмкость барьеров, г - представляет собой со- 
противление цепи, С„ - емкость воздушной среды 
с вкрапленными каплями воды, Я - сопротивле- 
ние, характеризующее ток проводимости при на- 
личии в разрядном промежутке субстанции, кото- 
рая состоит из воздуха и капель воды. 

Разобьем алгоритм определения параметров 
двухполюсника - эквивалентной электротехниче- 
ской модели озонатора в отсутствии разряда на три 
этапа. 

Первый этап - определение эквивалентной ём- 
кости С э и г на основе осциллограмм напряжения и 
тока при отсутствии капель воды в разрядном про- 
межутке. На этом этапе моделирования можно 
считать, что воздушная среда является идеальным 
диэлектриком и в цепи присутствует эквивалент- 


ная ёмкость С э 


с б с вв 

С б +С вв ’ 


образованная после- 


довательным соединением двух ёмкостей: С 6 - ём- 
кости барьера и ёмкости межэлектродного воздуш- 
ной промежутка С„, и сопротивление цепи г. 


Второй этап - проводится для разделения ём- 
костей С 6 и С„, которые определяются на основе 
расчёта распределения потенциала и напряжённо- 
сти электростатического поля электродной систе- 
мы с барьером. 


После определения величин г, С 6 , С„ наступает 
третий этап для расчета К и С т . Появление капель 
воды в разрядном промежутке вызывает появление 
проводимости §= \/Я, а также изменение емкости 
разрядного промежутка С вв . Для определения этих 
величин рассчитываются энергетические оценки с 
использованием осциллограмм тока и напряжения 
в случае, когда в воздушном межэлектродном про- 
странстве присутствуют капли воды. 


1. Постановка задачи для первого этапа 

Если действующее напряжение в цепи имеет 
импульсную форму, то переходный процесс, проис- 
ходящий в разрядном промежутке, удобно описы- 
вать с помощью интеграла Дюамеля [4, 5]. Приве- 


дем пример записи интеграла Дюамеля для опреде- 
ления тока і(і) через переходную проводимость §(1) 
и напряжение и(і), когда действующее напряжение 
и{і) определено на некотором временном интерва- 
ле, то есть и(1 ) является финитной функцией: 


К 0 = «(0) 8(1) + \ С ^ Ы у^-8(1 - т)с1т. 


СІТ 


( 1 ) 


Интеграл Дюамеля (1) - это решение прямой 
задачи определения тока і(і) по известным величи- 
нам и(1) и §(Г). В нашем случае представляет инте- 
рес обратная задача, задача определения §(1) пере- 
ходной проводимости по измеренным данным і(і) 
и и(1). Определение проводимости^/) эквивалент- 
но определению параметров схемы, так как §({) яв- 
ляется их однозначной функцией. 

Уравнение (1) это интегральное уравнение 
Фредгольма второго рода относительно переход- 
ной проводимости §(1). 

При отсутствии разряда и капель воды в разряд- 
ном промежутке в схеме замещения участвуют ём- 
кость барьеров С б , ёмкость воздушной среды С„ и 
сопротивление схемы г. В соответствии с приве- 
дённой схемой замещения (рис. 1 , б) переходная 
проводимость в интеграле Дюамеля равна: 

8і(0 = -е р1 , ( 2 ) 

г 


1 

где р = -—. 

С э г 

То, для определения параметров цепи г, С м , и С 6 
необходимо решить интегральное уравнение (1). 


Модельный пример определения параметров 
электротехнической модели озонатора 


Приведём пример расчёта параметров схемы 
электротехнической модели озонатора на основе 
интегрального уравнения Фредгольма второго рода 
(1). Для входного напряжения возьмём функцию, 
имеющую вид типичных осциллограмм напряже- 
ния при барьерном разряде (рис. 2, й): 


и(і) = 


1,5-10 


4 . е -« 5 У-у) в 


1 + р-( 1 -у) 


( 3 ) 



а 

Рис. 1. Эквивалентная схема озонатора в отсутствии разряда 
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где «=2,8- ІО 6 ; /5=8- ІО 5 ; 7 = 7/2 - коэффициенты не- 
линейной регрессии, определённые для реальной 
осциллограммы напряжения, Т — длительность 
импульса входного напряжения, здесь 7=1 мкс. 

Выберем в качестве модельных величин цепи сле- 
дующие параметры С г =5- 10 11 Ф, С и =3,372-10 -10 Ф, 
г=5 Ом. В соответствии с заданными величинами, ис- 
пользуя ур. ( 1 , 2 ), находим переходную функцию: 

§(!) = 0,2-ехр(-4,593-10 9 7). 



а 



б 

Рис. 2. Осциллограмма импульса воздействующего напря- 
жения и входного тока 

Находим ток в цепи, используя интеграл Дюа- 
меля (1). Результат вычислений, с использованием 
квадратурной формулы прямоугольников для вы- 
числения интегралов и кубических сплайнов для 
вычисления производных напряжения, приведён 
на рис. 2, 6. Теперь решаем обратную задачу - зада- 
чу определения переходной функции из интеграль- 
ного уравнения Фредгольма (1) с последующим 
вычислением всех электротехнических параметров 
цепи. Алгебраизуем ур. (1) 

N 

'к = и 08 к +Аі ^ ^ и к,т8щ > ( 4 ) 

/ 72=0 

здесь сіи к т І т) . 4 = /( 4 ); <? к =§(і к ); щ=Щ 0 ); 

с/1 

А(=Т/И— шаг дискретизации; т, кеО...М, /V — чи- 
сло точек дискретизации временного интервала. 

Запишем решение ур. (4) в матричной форме, 
введя обозначения, гіі ]={с!и кт }, Е - единичная ма- 
трица, С={&} - матрица проводимости, І={ 4 } - ма- 
трица тока: 


К = (ГО+^-Е, С = К _1 I. (5) 

А/ 

При обработке натурных экспериментов изме- 
ряемые величины имеют погрешности, поэтому к 
функциям и(!) и і(і) добавлен шум величиной 5 % 
от их максимального значения. На рис. 4, а, приве- 
дён результат восстановления функции проводи- 
мости (ломанная кривая) и её модель (гладкая кри- 
вая). После операции сглаживания фильтром - 
скользящим средним, восстановленная кривая 
полностью совпала с модельной. Для определения 
параметров схемы прологарифмируем восстано- 
вленную сглаженную кривую проводимости §(/), 
рис. 3,б,з. затем, используя метод наименьших ква- 
дратов, определяем наилучшее линейное прибли- 
жение Л/)=а/+Ъ на интервале высокого отношения 
сигнал/шум, из которого определяем параметры 
цепи г и С э : 

С э = - — ; г = ехр(-б). 
га 




б 

Рис. 3. Модельная (гладкая кривая) и восстановленная ( из- 
ломанная кривая) функции проводимости 

Расчеты, проведённые по выше описанному ал- 
горитму, дали следующие результаты: г =4,936 Ом и 
С э =44,82 пФ. Ошибки полученных величин не 
превышают 3 % от истинных значений для сопро- 
тивления и ёмкости соответственно (/= 5 Ом, 

с с вв с б = 

ее ^ б 

2. Второй этап - расчёт электростатического поля 

В качестве исследуемой системы рассмотрим 
электродную конфигурацию натурной лаборатор- 
ной установки, которая используется для очистки 
воды - группу круглых, равноудалённых металли- 
ческих цилиндров диаметром 5 мм, расположен- 
ных параллельно (рис. 5, размеры указаны в мм). 
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а б 

Рис. 4. Плоская электродная система: а) без барьеров ; б) с барьерами 


Металлические цилиндры вставлены в плотно 
прилегающие стеклянные колбы толщиной 1 мм 
(барьеры из кварцевого стекла, е=4, см. рис. 4, б). 
На каждый цилиндр подаётся напряжение 
Н=10...15 кВ с чередующимся знаком, при этом 
разность потенциалов на электродах составляет 
Н=20...30 кВ. Напряжение на электроды подаётся в 
виде импульсов длительностью 0,6 мкс. В верти- 
кальном направлении, через электродную систему, 
пропускаются капли воды различной фракции, не 
превышающие размер воздушного промежутка. 

Приведём расчет поля электродной системы 
при отсутствии диэлектрических барьеров 
(рис. 4, а). Для расчета поля такую электродную си- 
стему можно считать бесконечной в продольном 
направлении (в приближении электродов беско- 
нечной длины) и расчет осуществлять только для 
одной её секции, под которой будем понимать 
область, ограниченную плоскостями А и В. Напря- 
жение, приложенное к электродам, здесь и далее 
принимается 11=20 кВ. 


Определим функцию Грина - поле, порождён- 
ное точечным зарядом, находящемся между двумя 
заземлёнными плоскостями [6, 1\. 


Щх,У ) 


1 

2 Я8 0 



„;_,*(* + ■ > У + Уо )1 

1п 

21 

,^ 2 (х-х 0 ,у-у 0 ) 


2 1 1 


( 6 ) 


здесь і{х,у)=х+]у, и ;>о,, у 0 - координаты расположе- 
ния заряда, г 0 =8,85-1(Г 12 Ф/м. 


Для определения поля цилиндра между зазе- 
млёнными плоскостями необходимо решить инте- 
гральное уравнение Фредгольма первого рода, 
ядром которого является выражение (6). 


Ѵо = | (7) 


сіі 


где о(х,у) искомая плотность поверхностных заря- 


дов на цилиндре, 


сіз(і) 

сіі 


(іх(і) 2 сіу(і ) 2 

+ =К Ц - 


сіі 


сіі 


радиус цилиндра, ІІ 0 - потенциал цилиндра. 

Осуществляя алгебраизацию интегрального 
уравнения (7), получаем систему линейных алге- 
браических уравнений относительно ст 


<т=А- 1 ІІ 0 , (8) 

здесь А , у =Щх(1),у ( /,) ) /?„Д /; А і— шаг дискретизации; 
/,/е 0...А, N - число точек дискретизации контура. 
Выражение (8) является некорректным линейным 
алгебраическим уравнением [8], для решения кото- 
рого необходимо осуществить регуляризацию: 

ст=(А*А + Я-Е)“ 1 А*Н 0 , 


где Е - единичная матрица, А* - эрмитово-сопря- 
жённый оператор, Я=10 6 - параметр регуляриза- 
ции [8]. После определения а(х,у) заряда цилиндра 
определяем поле в электродном промежутке: 
(р(х 0 ,У 0 ) = 


Х.Л 

I а{х{1),у(1))\п 


1 


Ф(()-х 0 ,у(()-Уо). 


СІх(і) 


сіі 


сіі. (9) 


Осуществляя дискретизацию интегрального 
выражения (9) методом прямоугольников, опреде- 
ляем распределение потенциала Щх,у)=Яе(ср(х,у)) 
и функции потока Ѵ{х,у)=\т{ср{х,у)) в дискретном 
множестве точек. Подействовав оператором V в 
комплексной форме на выражение (9), можно 
определить напряжённость поля в любой точке 

Е(х,у) = —^(р(х,у) = 

02 

_ 2 С а(х(1),у(()) 

о 2 (х(1)-х 0 ,у(1)-у 0 ) сіі ’ (10) 


здесь черта над і - знак комплексного сопряжения. 

Теперь можно определить поле электродной си- 
стемы при наличии барьера (рис. 5, б). Разделив 
полученный вектор Е(х,у) на е { в области простран- 
ства, заполненного воздухом, и на е ъ где простран- 
ство заполнено стеклом, получим СЕ х (х,у) и 
СЕ 2 (х,у) напряжённости в соответствующих обла- 
стях с точностью до мультипликативной постоян- 
ной - С. Для определения постоянной С запишем 
выражение для разности потенциалов І7 0 через на- 
пряжённость поля с учётом (10): 


| Е х і(х,0)с1х + 1 Е х2 (х,0)сіх 

ѵ*0 *1 у 
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в) Линии равного потенциала ІІ(х,у) = сопеі и равного потока Ѵ(х,у) = сопзі 



Рис. 5. Картина распределения потенциала ІІ(х,у) и напряжённости электрического поля Е(х,у) в электродной системы с барье- 
ром. Для сравнения пунктиром нарисованы профили значений напряжённости и потенциалов поля без барьеров 


здесь Е ХІ (х, 0) и Е Л (х, 0) х-овые компоненты вектора 
электрической напряжённости в первой и во вто- 
рой средах соответственно, и у 0 - точки простран- 
ства с потенциалом, равным нулю, х х и у х — точки 
пространства с потенциалом, равным II п . После 
определения константы определяем потенциал в 
любой точке межэлектродного промежутка ср(х,у). 
х У 

<р(х,у)=- | Е х (х',у)сіх'- | Е у {х,у')ііу'. 

хо уо 

Результат восстановления представлен на 
рис. 5. Там же представлены профили модуля элек- 
трического поля при х=0 и у= 0. 

Расчет ёмкости электродной системы 

Геометрия электродов рассматриваемой систе- 
мы с барьером имеет неканоническую форму, а 
межэлектродное пространство состоит из двух ку- 
сочно-однородных сред: воздушной - е,=1 (при 


отсутствии капель воды) и заполненной кварцевым 
стеклом г 2 =4. Расчёт ёмкости такой системы мож- 
но осуществить с помощью рассчитанного элек- 
трического поля. Технология расчета такова: для 
рассчитанного электростатического поля исполь- 
зуется теорема Гаусса - определяется заряд какого- 
либо одноимённого электрода в одной секции. 


2 п 

с 1 = е о\ г(х(1),у(1))[Е х (х(1),у(1)) 

О 


< 1 у( 1 ) 

сіі 


-Е (х(1),у(1))^^-\сІІ, 
аі 


здесь Е х (х,у) и Е у (х,у) - компоненты вектора элек- 
трической напряжённости, х(1) и у(1) - параметри- 
ческое уравнение контура, обхватывающего элек- 
трод (возможные контуры приведены на рис. 6 
пунктиром). При известной разности потенциалов 
ІІ а =(р-(р 1 ёмкость электродной системы определя- 
ется по формуле [9]: 
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С = 


Ц- 


(П) 


Выражение (11) использовалось для каждой 
области, занимаемой стеклом и воздухом в отдель- 
ности. Разность потенциалов для воздушной обла- 
сти определялась как 11 а =(р-(р ъ где ф,=0,87- ІО 4 В - 
потенциал границы раздела сред (рис. 5,6- про- 
филь распределения потенциала). Для области с 
барьером разность потенциалов определялась как 
11 п =(р { -(р2 и составляла величину 6^2=0, 1 3- ІО 4 В, где 
<р,= 1Т0 4 В и <^=0,87-1 0 4 В. В результате были полу- 
чены значения ёмкостей на единицу длины для 
воздушной и заполненной стеклом (барьером) сред 
соответственно: 

С в =5, 039-КГ 11 Ф/м; С б = 3,372-1 0' 10 Ф/м. 

Ёмкость всей секции определяется как после- 
довательное соединение ёмкостей сред: 

С 4 384- ! (Г 11 Ф/м. 

с б + с в 


Учитывая, что мы имеем четыре секции, и ём- 
кости секций соединены параллельно (рис. 6), ём- 
кость всей системы определяется как параллельное 
соединение ёмкостей секций: 

С« ч =Сз+С э +С э +С 3 =1,74-10- 1 “Ф/м. 



Рис. 6. Эквивалентная электрическая схема секции разрядника 


напряжения и тока, определим напряжение на ем- 
кости водо-воздушной среды С„„ по второму закону 
Кирхгофа 


«Сев (0 = «(О 


/ 

і{і)г + 

V 


'-'б у 


Рассмотрим временной интервал, равный удво- 
енному значению длительности импульса входного 
напряжения (рис. 7). По окончании переходного 
процесса суммарная энергия элементов схемы яв- 
ляется чисто активной и составляет сумму энергий, 
потребляемых на сопротивлениях г и К. Вклад ре- 
активных составляющих уменьшается и становит- 
ся незначительным (пунктирная и штрихпунктир- 
ная линии на рис. 7), поэтому для конечного вре- 
менного участка можно записать 
]Ѵъ]У г + ]Ѵ к . 



Рис. 7. Распределение суммарной энергии, активных и реак- 
тивных её составляющих: 1Л/„ И/, - активные энергии, 
потребляемые на сопротивлениях г и Р; I Ѵ Сб , I Ѵ Сш - ре- 
активные энергии; И/ - полная энергия, потребляе- 
мая схемой 


3. Определение параметров схемы замещения 
при наличии капель воды 

При наличии капель воды в электродном проме- 
жутке определить емкость системы становится за- 
труднительным. Единственным способом в этом 
случае является метод описания электрического раз- 
ряда как объекта электрической цепи [1-4], исполь- 
зуя осциллограммы входного тока и напряжения. 

На первом и втором этапах были определены 
величины, являющиеся постоянными С 6 и г, най- 
денные по вышеизложенным алгоритмам на осно- 
ве расчёта поля электродной системы и осцилло- 
граммам напряжения и тока при отсутствии разря- 
да и капель воды в разрядном промежутке. Приве- 
дём алгоритм определения изменяющихся параме- 
тров цепи К и С„, величина которых зависят от ко- 
личества, фракции и загрязнённости капель воды, 
находящихся в межэлектродном промежутке, ра- 
спределения капель по объёму разрядного проме- 
жутка (рис. 1, а). Зная осциллограммы входного 


Так как найденное напряжение равно напряже- 
нию на искомом сопротивлении К и к (Г)=и Сіі (1), 
рис. 1, а, то энергия на сопротивлении определяет- 
ся выражением: 

і , і 

\Ѵ к =\Ѵ-\Ѵ г = к\і к 2 (і)6і = -\ и к 2 ( 13 ) 

• К і 

о о 


Откуда Я = 


1 


IV - IV г 


]и к 2 №. 


а емкость водо-воздушной среды определим из из- 
вестного соотношения 


г '(0 = с вв 


бы к (() 
бі 


с вв = 1 ( 0 - 


6и к ( I ) 
бі 


-1 


. (14) 


Пример решения модельной задачи 

Рассмотрим модельную задачу восстановления 
К и С ав . В качестве входного напряжения выберем 
импульсную функцию, заданную выражением (3). 
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Выберем в качестве модельных величин цепи, изо- 
бражённой на рис. 1, а, следующие параметры 
С,=5-10- п Ф, С= 337,2 пФ, г=5 Ом, Л=1000 Ом. В 
соответствии с заданными величинами, решаем 
прямую задачу - определяем переходную функцию 
§ 2 (1) для входного тока цепи второго порядка: 

82 (0 = с б ІАРі е х Р Ш) + А іРі ехр (р 2 *)} = 

= А схр( р\і)+ В 2 ехр (р 2 1). 


где А І и А 2 определяются из системы уравнений: 

Р2 гС 6 + 1 . 


А А 2 + 1 - 0 ; 

1 

СбіАРі + Ар 2 ) = — ; 

г 


А = ~ 


(р 2 -рі)гС б 


А аДкіі , (і5) 

( Р2~Р\У С б 


а характеристические корни р { и р г , входящие в си- 
стему (15), определяются из квадратного уравнения: 

С б С вв Кг р 2 +(С вв К +С б г +С б К)р + 1=0. 

Теперь, при известных входном напряжении - 
и(1) и переходной функции § 2 (1), используя инте- 
грал Дюамеля, получаем «осциллограмму» входно- 
го тока. Энергия в цепи распределяется в соответ- 
ствии с выражением (13), из которого находим со- 
противление К, а параллельное соединение эл- 


ементов К и С„ позволяет определить величину ём- 
кости С„„, используя соотношения (14). 

По вышеописанному алгоритму определено, 
что Д=998 Ом, С е = 338,7 пФ, что не превышает 3 % 
ошибки от заданных значений. 

Заключение 

В работе предложен алгоритм синтеза электри- 
ческой цепи - электротехнической схемы замеще- 
ния барьерного разряда. Алгоритм основан на рас- 
чёте распределения электрического поля с после- 
дующим вычислением ёмкостей электродной си- 
стемы с барьером в случаях, когда в межэлектрод- 
ном воздушном промежутке отсутствуют капли во- 
ды. При наличии воды в воздушном межэлектрод- 
ном пространстве алгоритм расчета параметров 
электрической цепи, основывается на решении 
интегрального уравнения Фредгольма, связываю- 
щего ток и напряжение с переходной функцией 
проводимости цепи. Искомая проводимость содер- 
жит в себе все параметры электротехнической мо- 
дели озонатора. На примере численного экспери- 
мента продемонстрирована удовлетворительная 
работа алгоритма по восстановлению параметров 
электротехнической схемы замещения озонатора. 
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